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Abstract 

We consider a quantum two-dimensional O(N) (8> 0(2)/0(N — 2) (g> 0(2)di a g nonlin- 
ear sigma model for frustrated spin systems and formulate its 1 /A-expansion which 
involves fluctuating scalar and vector fields describing kinematic and dynamic in- 
teractions, respectively. The ground state phase diagram of this model is obtained 
within the l/TV-expansion and 2+e renormalization group approaches. The temper- 
ature dependence of correlation length in the renormalized classical and quantum 
critical regimes is discussed. In the region p- m < p ut 5 Xin < Xout of the symmetry 
broken ground state (/9i n ,out and Xin.out are the in- and out-of-plane spin stiffnesses 
and susceptibilities), where the mass M u of the vector field can be arbitrarily small, 
physical properties at finite temperatures are universal functions of Pxnaat: Xin.out > 
and temperature T. For small these properties show a crossover from low- to 
high temperature regime at T ~ M u . For p m > p out or Xin > Xout finite-temperature 
properties are universal functions only at sufficiently large Ma- The high-energy be- 
haviour in the latter regime is similar to the Landau-pole dependence of the physical 
charge e on the momentum scale in quantum electrodynamics, with mass play- 
ing a role of e" 1 . The application of the results obtained to the triangular-lattice 
Heisenberg antiferromagnet is considered. 
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1 Introduction 



The description of frustrated systems is a long-discussed subject [Tf2"f5] . These 
systems pose an important and interesting problem of condensed matter physics, 
and also provide a test of various methods of quantum field theory [HI51I6] . 
Competing interactions leading to frustration favor strong quantum fluctua- 
tions and therefore tune the system towards quantum phase transition (QPT) 
into the phase without spontaneous symmetry breaking. 

Examples of frustrated systems are the Heisenberg triangular lattice antiferro- 
magnets (TLAF). In particular, VCI2 and VBr2 are layered compounds with 
the triangular lattice structure within a layer. The corresponding Heisenberg 
model with nearest-neighbour exchange interaction is shown to have a non- 
collinear ordered ground state with the sublattice magnetization suppressed 
by quantum fluctuations [7.8 9J. TLAF can be further tuned to QPT into the 
spin-liquid state by including additional next-nearest [10J or ring exchange in- 
teractions [UJ, impurities or charge fluctuations [T2]. Another class of systems 
which show a non-colinear (helimagnetic) order are rare earth elements [13J. 

Previous investigations of frustrated and non-collinear antiferromagnets mainly 
concentrated on the critical behaviour near the magnetic phase transition. 
The nonlinear sigma (NLcr) model [Sj and Landau- Ginzburg- Wilson (LGW) 
approach [HfTi] were applied to this problem. These approaches predict criti- 
cal exponents which are different from those of the standard 0(3) universality 
class. Recent investigations within the non-perturbative renormalization group 
approach predict, however, a first-order phase transition in three dimensions 
[6]. The latter method was also applied to the description of temperature 
properties near the quantum first-order transition in two-dimensional (2D) 
frustrated antiferromagnets [TB] . 

The consideration of physical systems requires a description of a broad temperature- 
and control parameter range. Even not too close to QPT one can differentiate 
the high- and low- energy degrees of freedom. While the latter correspond to 
"slow" degrees of freedom and can be described by a continuum model, the 
former, "fast" degrees of freedom, can be absorbed into renormalization of 
parameters of the model. Thermodynamic properties at low temperatures can 
be expressed in this way as universal functions of ground state parameters. 
The corresponding functions describing quantum antiferromagnets can be ob- 
tained using the above mentioned NLcr model. In the classical case this model 
corresponds to a continuum limit of the Heisenberg model. For the collinear 
quantum antiferromagnets the NLcr model was first derived by Haldane in the 
framework of the l/S'-expansion p^6|TT7] . Despite the way of derivation, this 
model is applicable in both the symmetry broken and the symmetric phases. 
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The finite temperature properties may become non-universal in certain cases, 
e.g., near first order phase transitions, and also for systems described by the 
quantum field theories at and above their upper critical dimension. The latter 
theories are non-renormalizable and therefore lead to non-universal properties. 
Especially interesting example of such a possibility occurs when the theory 
contains two types of interaction terms, one of which is below and another is 
above their upper critical dimension. While the former interactions are infrared 
relevant and produce universal contributions to physical properties, the latter 
produce non-critical but non-universal contributions. As it is shown in this 
paper, such a situation is realized in frustrated spin systems, in particular 
TLAF. 

The noncollinear magnetic order in frustrated systems is described by so-called 
0(3) <S> 0(2)/ 0(2) NLa model [IB]. This model includes only terms appear- 
ing in the large-S* limit and is therefore quasiclassical. In comparison to the 
0(3)/ 0(2) model for the square lattice this model has different 0(3) <8> 0(2) 
symmetry and a matrix order parameter [5]. For the description of physical 
systems the above discussed quasiclassical model should be extended to include 
all the relevant terms which turn out to be responsible for non-universality in 
certain parameter range. The construction of the quantum 0(3) <S> 0(2)/0(2) 
NLcr model allows one to investigate temperature behaviour of physical prop- 
erties (e.g., correlation length) above the symmetry broken ground state and 
near QPT. The implementation of this task, however, calls for the application 
of non-perturbative methods to deal with the proximity to symmetric phases 
[19] and critical fluctuations [2"5] . 

An useful non-perturbative tool of investigating the low-temperature proper- 
ties of NLcr models of spin systems is the 1/iV-expansion, N being the number 
of the spin components (N = 3 for the physical case) [2U[21|l22ll23ll24|l25j . 
Contrary to the field-theoretical renormalization group (RG) approach, this 
method does not suppose universality and renormalizability of the model. For 
classical case this method introduces a single constraint J2i n ? — A/" (n is the 
fluctuating unit length order parameter, M is the number of sites) which is 
exact in the limit N — > oo [26f27j . In spite of fluctuations of spin length, 
this approximation appeared to be surprisingly good for description of ther- 
modynamic properties of collinear antiferromagnets [22,28,29,30]. Considering 
1 /^-corrections allows to obtain in addition a correct description of the criti- 
cal region. Note that previous considerations of the thermodynamic properties 
of frustrated antiferromagnets within the 1/iV-expansion in the CP^" 1 formu- 
lation were constrained to a state with a higher 0(3) x 0(3) symmetry [2J. 

The 1/iV-expansion is to some extent complementary to the RG approach. In- 
deed, the 2+e RG method is perturbative in coupling constant and constrained 
to the vicinity of two space-time dimensions, but it is non-perturbative in the 
number of the components. In contrast, the 1/iV-expansion is non-perturbative 
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in the coupling constants and has no restriction on the dimension. Therefore we 
consider both the methods in our study of frustrated spin systems to have an 
occasion of comparing their results. Simultaneous using these methods enables 
one to determine physical properties in a broad temperature and parameter 
region. 

The content of the paper is as follows. In Sect. 2 we discuss the quasiclassical 
0(3) ®0(2)/0(2) NLcr model and its extension to the quantum case. In Sect. 
3 the generalization of the quantum NLa model to the iV-component order 
parameter and a possibility of performing its 1/iV-expansion are considered. 
We examine the saddle point N = oo approximation and consider its stability 
in Sect. 4, then obtain fluctuation corrections to the Green's functions in Sect. 
5. In Sect. 6 we study the renormalizability of the model in the ground state 
and obtain the RG-flow phase diagram in both the 1/iV-expansion and 2 + e 
RG approaches. In Sect. 7 we apply the 1/iV-expansion (and the RG approach 
where it is possible) to calculate correlation length in the renormalized clas- 
sical and quantum critical regimes. The application of the results to TLAF is 
considered in Sect. 8. Our results are summarized and discussed in Sect. 9. 



2 Effective action 

As a prototype of frustrated systems, we consider TLAF described by the 
Heisenberg Hamiltonian 

Hfteis = J X! S?.Sj, (1) 

(iJ) 

where J > is the antiferromagnetic exchange parameter, S* are spin oper- 
ators on the sites of the triangular lattice, and denotes the summation 
over nearest neighbours. Even in the classical case, due to frustration, there 
are no configurations that would minimize energy of exchange interactions 
between all nearest neighbours. This results in the noncollinear ground states 

(Si) = S(e 1 cos(Qxi) + e 2 sin(Q Xi )), (2) 

where Q = (27r/3, — 27r/v^3) is the wave vector of the AFM structure, X j is the 
radius vector of z-th site, and ei, e 2 are two orthonormal vectors (ef = e| = 1, 
ei • e 2 = 0) which determine the plane of spin alignment. An example of such 
structure is shown in Fig. 1. 

Quantum fluctuations lead to the suppression of the sublattice magnetization 
(Sj)o — > Z(Sj)o (Z < 1, Z = 1 in the classical limit). For example, in the 
extremely quantum case 5 = 1/2 the result of SWT in the second order in 
1/S is Z = 0.53 [31] (for non-frustrated square lattice with S = 1/2 the 
corresponding result is Z = 0.61 [32] ) 
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Fig. 1. The classical ground state magnetic and corresponding plaquette long-range 
orders 



The effective low-energy and long- wave-length field theory for the Hamiltonian 
(CQ) was derived in Ref. [18] in the quasiclassical limit S —>■ oo using the Haldane 
mapping [T61ITT) . The corresponding action has the form of a NLcx model 



Sbr = 7^— f 1/T dr ( d 2 x [|^ ei | 2 + \d T e 2 \ 2 + |<9 r e 3 | 2 + c 2 n (|V ei | 2 + |Ve 2 p 
2c m g Jo J 1 

(3) 

where g = 6a/ S is the coupling constant (a is lattice parameter), orthogonal 
vectors ei, e 2 now acquire space-time coordinate dependence, e 3 = [ei x 
e 2 ], Qn is the in-plane spin-wave velocity. It is assumed that there is some 
regularization procedure for the wave vectors and frequencies (e.g. using a 
cutoff parameter A). 

Since quantum effects are not small for physical systems, other terms in the 
action ([3]) may be present, besides the renormalization of the spin- wave veloc- 
ity Qn and the coupling constant g. One can obtain corresponding quantum 
action using symmetry properties of the original Hamiltonian (JTJ) with respect 
to the SO (3) rotations in the spin space, translations and the time reversal 
r — > — r, ei — > — ei, e 2 — ► — e 2 , as well as the transformation 



ei — ► ei cos ip + e 2 sin (p 
e 2 — >• — ei sin (p + e 2 cos <p, 



where ip = 2irm/3 (m is an integer number); the latter symmetry follows from 
the invariance of the model ([1]) with respect to the translation by the lattice 
vector. The action, which is invariant under these operations, has the form: 
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1 fl/T 
s = - 

2 Jo 



J drj d 2 x[p 1 (\d T e 1 \ 2 + \d T e 2 \ 2 )+p 2 \d T e 3 \ 2 + (5) 
+ gi(|Vei| 2 + |Ve 2 | 2 ) + g 2 |Ve 3 



where we have included only the second-order derivative terms. Coefficients 
Pi >2 and gi j2 are related to bare spin stiffnesses and susceptibilities of the model 
§Q (see Sect. 3). The action ([5} neglects the contribution of the Berry phase 
term, which is expected to be not important for the ordered ground state and 
low T [33]. Indeed, singular configurations e a are separated by an energy gap 
of order of J from the ground state. 

The effective action §5§ can be related to the so-called massive CP*'" 1 model 
proposed in another description of frustrated systems in terms of the spinon 
field, connected to the order parameter field e a by the local two-valued func- 
tion [2]. However, due to the two- valued character of this mapping, the massive 
QpM-i moc | e i h as t be accomplished by the Z 2 lattice gauge field, which pro- 
vides confinement of spinons [31]. In the following we will consider the field 
theory (jSJ) which describes spinons in the confined phase and do not consider 
the possibility of their deconfinement [35J. 



Generalization of the action to the iV-component order parame- 
ter 



To generalize Eq. (jS]) to the iV-component case we consider the action 



1 r l / T r r 

S g cn = -J dr J d 2 x\p 1 (\d T e 1 \ 2 + \d T e 2 \ 2 ) +p 2 (\d T e 3 \ 2 ^ h |^ejv| 2 ) 

+gi(|V ei | 2 + |Ve 2 | 2 ) + g 2 (|Ve 3 | 2 + • ■ ■ + |Ve^| 2 )] , (6) 

with e a (a = 1, . . . , N) being iV-component vectors forming the orthogonal 
basis. We follow the construction of the 1/iV-expansion in the collinear case 
[25ll28ll33j . i.e., introduce Lagrange multipliers X a/ 3 for the constraints e^e^ = 
5 a/ 3 and calculate the Gaussian integral with respect to e a . Since, contrary 
to the collinear case, \ a p is a iV x iV-matrix, even the N — > oo limit is non- 
trivial: it is necessary to sum up an infinite subsequence of the so-called planar 
diagrams, which is a rather difficult task [2~TH2"o] . 

To formulate an appropriate 1 /iV-expansion we consider first the case q 2 = 0, 
p = p\ = p 2 which holds for large enough spin S. In this case it is possible 
to integrate out the vector fields e 3 , . . . e^ exactly (note that only ei and e 2 
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appear in Eq. (T5D) and obtain the effective action (see Appendix A): 



1 fl/T 

'S'qeaf = ^ J Q ^ T 



J d 2 x [Jf (|^ ei | 2 + |^e 2 | 2 ) - 2J£|(e 1 , d M e 2 )| 2 ] , (7) 



where /i = r,x,y, J\ = 2p, Jf = Jf = q 1 , J 2 = p, and Jf = Jf = 0- 
The action (j7j) is justified also for N = 3, which can be verified by direct 
substitution e 3 = [ei x e 2 ] into Eq. (j3J); its classical analog was studied in Ref. 
[5]. Since ([7]) contains only two fields, the 1/iV-expansion can be constructed 
for this action (see SectSJ). Note that in both the cases Jf < and J 2 > 
the action ([7]) is stable (in the latter case due to the presence of the nonlinear 
constraints). 



The above-discussed procedure of integration over the vector fields e 3 , . . . , 
cannot be easily generalized to arbitrary pi j2 and gi i2 . In this case we find that 
the effective action contains not only terms with squares of gradients (as in Eq. 
(J7|)), but also additional terms with higher gradient powers. We assume that 
these terms give subleading contribution to physical properties and restrict 
ourselves to the action (J7|) with J* = (a = 1, 2) and arbitrary J£ (/x = x, r), 
referred in the following as quantum effective action for frustrated systems 
(QEAF). This action is a generalization of the quasiclassical action considered 
by Azaria et al. p5j. The coupling constants are related to the bare spin 
stiffnesses p? n out and susceptibilities Xi n out 



V — v° T x — n° 



out ' 



Jl = Xl - V2X° n , J2 = Pout " V2p° n 



where the subscripts 'in' and 'out' correspond to the in-plane and out-of- 
plane modes [31][2]. QEAF has a form of 0{N) ® 0(2)/0(N - 2) <g> 0(2) diag ~ 
0(N)/0(N — 2) NLcr model, 0(2)-symmetry being defined by the relations 
(jl]) with arbitrary (p. Hence, we find two types of terms in QEAF. The first 
one has the same structure as in the original action (jH]). The second (current) 
term introduces a dynamical quartic interaction additionally to the kinematic 
one, the latter originating from the nonlinear constraints. This term is related 
to the square of conserved Noether' current J7 12 = (J^ 1 — Jf )(ei, <9 M e 2 ) for 
0(2)-symmetry. 
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4 1/iV-expansion and solution at N — > oo 



To construct the 1/iV-expansion, it is convenient to decouple quartic term 
(e 1 ,<9 At e 2 ) 2 introducing the Hubbard- Stratonovich field A^. 

£ = ye Q X a %-^S P A, (9) 
X aP = -( Jf - JZ)dl5 aP + JZ[(td, - a^A^Xp + A Q/3 , (10) 

where cr^g is the Pauli matrix, 

and An,Ai2,A22 are the Lagrange multipliers for the constraints ei • e 2 = 
0, ei 2 = e 2 2 = 1. Integration over A^ is performed along the real axis for 
J 2 > and along the imaginary axis for J 2 < 0; the integration over A is 
performed along the imaginary axis. To deal with the symmetry broken phase 
it is convenient to represent the field e a as a sum of a fluctuating field e a and 
the uniform static order parameter m a which is chosen to satisfy / dxe a = 
m Q / dx. Integrating out the field e a we obtain the following representation 
for the partition function: 

Z = J D[X af3 ] J D[Afj] J D[m Q ]exp(-5 cff [A,A,m,I,h]), (12) 



N 

5 cff [A,v4,m,j,h] = — splnX (13) 



— / dr I d 2 x 



2 

Sp X(x) + m a ip a 0mp + — h a m 
2, Jo J L iv 

—^da/N - (fasms^X'Xpijp/N - ^/3 7 m 7 ), 

where we have introduced the sources j Q and h a of the fields e a and m Q , 
respectively, and 

ip af3 = J^A^Sap - Jgialpd^A^} + \ aP . (14) 

' Sp' denotes the trace over the index a = 1,2, ' sp' includes also an integral 
over the space-time variables. 

Calculating the derivatives of Eq. ( fT2l) over the source fields we obtain the 
general expression for the Green's function G^p = (^e l a (x)e^(y)\: 

G% = (jjiX- 1 )^ + [(rVm% ® (X- W)/j] + mX) > (15) 
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where <g> stands for the tensor product with respect to coordinates, the av- 
erage is taken with the weight exp (— S e ff). From the field-theoretical point 
of view, this Green's function corresponds to the propagation of elementary 
excitations of the symmetry group O(N) x 0(2) with quantum numbers 1=1, 
t = ±1 corresponding to the 0(N) and 0(2) subgroups, provided that these 
symmetries are not broken. 



In the N — > oo limit the path integral ffl2|) is dominated by the vicinity of 
saddle points of the functional S e s[\, A,m], if they exist. The saddle-point 
equations have the form 

A A A 

Scs = 0, ——S eff =0, —S eS = 0. (16) 



5\{x) e " ' 5A{x) vu ' 6m 

The coordinate independent saddle point is \ a p = v 2 5 a/ 3, A T = A x = A y = 0, 
mi • m 2 = 0, mf = m 2 , = (Tq, with do and v satisfying the equations 

d 2 k 



vtio = 



where 



o; n = 2imT are the bosonic Matsubara frequencies, and Tq = t; 2 / Jf . It is shown 
in Appendix B that ReS^ has its smallest value at the saddle point (TT7j) . as 
compared to other possible, in general space-time dependent, saddle points. 
Therefore the saddle point (IT7j) provides mam contribution to the partition 
function at iV — > oo. 

The equations (ITTj) coincide with those for the collinear case [28] and were 
shown to have a solution at T = with either o~ ^ (for large enough Ji), 
describing the ordered phase, or v ^ (for small enough J x ) correspond- 
ing to the symmetric (paramagnetic) phase. These two states are separated 
by a second-order quantum phase transition. In agreement with the Mermin- 
Wagner theorem only symmetric solution is possible at T > 0. The finite T 
region above the symmetry broken ground state is denoted as the renormal- 
ized classical regime, whereas the region above the quantum critical point is 
referred to as quantum critical regime [36J. 



5 Fluctuations 



To go beyond the saddle point approximation considered in Sect. HJ we need 
to take into account fluctuations of the A and A fields. Expanding the action 
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near the saddle point (v 2 , 0) up to the second order in fluctuations of these 
fields, we obtain: 



5 eff = j J ^n(g)(|«n| 2 + 2|a 12 | 2 +|a 22 | 2 ) 



N r d 3 a 

+ Y J J2^ D MMQ)M-Q), (19) 



where a = — i(X — v 2 ), 

n(g) = n(g) + 2a 2 G (q), 11(g) = J ^G (p)G (p - q) (20) 

and 

D/w(q) = D^uiq) - 2aoJ^J^q ll q u G (q) 

- - r d 3 n (21) 
D, v {q) = 2 J 2 M <W - « J J^3 G o(p) G o(p - 1)(2P» ~ qJ&Pv ~ to) 

are the inverse propagator of the fields A and A^. We use the notation p^ = 
(p, iu) n ) and / d 3 p/ (27r) 3 for TJ2 n I d 2 p/ (2ti) 2 (and similar for q). Hereafter we 
choose the units c 2 nl = 1 and use a regularization of the integrals u 2 + p 2 < A 2 . 
Note that, in contrast to the collinear case, the action ([7]) is not relativistic 
invariant, so that the spin- wave velocities will be renormalized by the 1/N- 
corrections. 

Using the expression for the effective action f[T§l) we obtain the Green's func- 
tion 



Ggj(p) 



SjjSg/3 1 

{ "" N Gu\p) + J 1( 5r2 + S(j9) (22) 
-C i $\p)-C$ ij (p) + 5( P )mimj 3 . 

The first term in Eq. (I2"2"j) is isotropic and non-zero even in the symmetric 
phase. The self-energy is given by 

E(p) = E 1 (p)+E 2 (p), (23) 

where the contributions 

v / x 3 /• rf 3 g G (p-g)-G (g) 

(27F 5S) (24) 

and 

E 2 (p) = ^ / -^T^J^quG^q) - (2p - q),(2p - q) v G,(p - q)]D~\q)^ 

(25) 
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come from the fluctuations of A aj g and fields, respectively. The second and 
third terms in Eq. (1221) are due to the spontaneous symmetry breaking, 

^S>) = GKp^ml (a aj (p)ap S (-p)) , (26) 
C ( S J (P) = ^ T Gl{p)J^r 2 p, Pu al^ 5 m} 1 mlb-\p)^ (27) 



These terms contribute to the transverse in-plane and longitudinal modes in 
the leading order in 1/N (see [29|). In particular, they provide the correct 
value for c; n , which is different from c out . The mass correction 5r and static 
order parameter m a will be calculated in Sect. 7 and Appendix [D] in the first 
order in 1/N. 

For T = we obtain the propagators of the A and A fields for g<C A 



TT , . 1 1 q 
n w) = 12- A — arctan—, 



D,M) = J -^f- \M^ V + L u - Ml ) [n(q) + 2J,al\ \ , (28) 



Ttiq) = arctan 

W 8vrg 2r 



where 



M» = -^(Ji - J$) + — 2 (29) 



is the mass of the field A^. The last term in Eq. ([29]) originates from our 
regularization scheme which does not preserve the gauge symmetry for J\ = 
J%. This term, however, describes the renormalization of coupling constants 
due to short wave-length fluctuations and therefore should be retained. For 
J 2 M > we have M M > provided that J x - J% > 0, i.e. xL Pin > °- The 
propagators at q < A (when they are needed) can be evaluated numerically. 

For J 2 < the field A^ = iB^ is purely imaginary and the inverse propagator 
of the field equals —D^ u (q). In this case we obtain M M < 0, except for 
extremely small J\ where Eq. (1291) is not applicable, since we have neglected 
r in comparison with A; such a small J\ are not considered in the following. 
For any finite A and J 2 < the propagator —D^^q), Eq. fT2~T]) . is positively 
defined (this property can be verified by the convolutions of —D^ u (q) with 5^ v 
and q^q u ), which is in agreement with the stability of the saddle point (IT7j) 
proven in Appendix B. In fact, the propagator —D^ u (q) is a monotonously 
decreasing function of q, having a minimum for q = A. In the isotropic case 
J 2 = J 2 the smallest eigenvalue ("energy gap") corresponding to q = A is 
2J^/| J 2 I + 0.015 163 A. Due to small numerical coefficient this "energy gap" 
can be small enough for 2J-f/|Jf | — > 0. 
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6 Renormalization in the symmetry broken state 



For practical use of formulae of Sect. O the dependence of the observable 
quantities on the ultraviolet cutoff A should be absorbed into redefined bare 
parameters of the continuum model. In the present formulation the proce- 
dure of renormalization expresses renormalized uniform susceptibilities Xm, out, 
stiffnesses Pi n , out, and the order parameter |(S)| = Zi|m a | {Z\ is some non- 
universal factor), a = 1,2 in terms of the bare parameters of QEAF. This 
renormalization procedure at T = is similar to that in the quantum field the- 
ory and removes also ultraviolet divergences at finite T [28129] . The infrared 
divergences in the vicinity of critical points can be removed by additional 
renormalizat ion . 

To perform renormalization, let us choose the vectors m; along the axes i 
(i = 1,2). Consider the Green's functions G^ N and Gff corresponding to the 
out- and in-plane modes, respectively. In the leading order in 1/N we have 
(see. Eq. fl22J)): 

= -Uo(p) = ^ (30) 

with Xout = Pout = GqJi "i Xout = Pout in our units. Substituting this into 

Eq. ([IT]) we obtain: 

+ < 31 > 

Another relation can be derived from the Green's function 

1 - 



hoi 



p/8 + 2o%h 



- 2 (M x + tt(p))lu + (M T + 7r(p))p 

-h<T - 



M T (M X + tx(p))uj 2 + M X (M T + ff (p))p 2 



(32) 



(see Eqs. ( 122]) . (1251) ). where tt(p) = p/16+2p out /N. For small momenta p <C p out 
we find the spin- wave pole of Gf^(p): 



GH(p^ Pout ) = ^ - 2 , (33) 

where 

y . = 2Xout D . = ^25* (M) 

Xm l + 2 X ont/(M-NY Pm l + 2p aut /(M*N) 1 ' 

are the in-plane susceptibility and spin stiffness. From Eq. (}3i]) we see that the 

masses M x = 2N~ 1 p out / (2p out /p in - 1) and M T = 2N- 1 x oa t/ (2Xout/Xin - 1) of 

the field are expressed through the observable quantities. Combining Eqs. 

291 ). ( 15T1) and ( 1341) we arrive at the result for 

J2 

J * = C,[J 1 -A/(2n 2 )]+A/(3n 2 y (35) 
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where C T = [1 - Xi n /( 2 Xout)] 1 and C x = [1 - p in /(2p out )] 1 . Supposing A- 
independence of Pi n ,out, Xm,out an d varying A we arrive at the renormalization 
group transformation of the action ([7]), determined in the limit N — > oo by 
the formulae (13"Tj) and (1351) . Note that renormalizations of J T a and are inde- 
pendent from each other. The generalization of the RG flow for arbitrary N 
is considered below. 

For < (i.e. J2 < 0) the denominator of Eq. (1351 vanishes for A = A£, 

A* = -^^=^(1^1-^) (36) 

(and similar for A^). In fact the RG flow which begins in the region J2 < 
cannot be extended to A > since — > —00 for A — > A^ — 0. This 
corresponds to inapplicability of the continuum model ([Zj) at J2 < and 
scales A > A(f, where one cannot find bare parameters corresponding to given 
renormalized Pi n ,out> Xin.out- This peculiarity of the flow arises due to increase 
of the mass |M M | by fluctuations (see Eq. ff29]) ). 

The restriction on the RG flow is similar to that in the 4 -model in D > 4 
and quantum electrodynamics where the values of the renormalized coupling 
constant (charge) g r > g c (A) are not physically accessible, <? C (A) vanishing in 
the limit of infinite cutoff [37]. For the frustrated systems the inapplicability 
of the continuum QEAF may correspond to nontrivial change of excitation 
spectrum similar to that earlier obtained for the TLAF within the series ex- 
pansion [38J and renormalized spin- wave theory [39]. Note that the scale 
can be small enough in comparison with the inverse lattice constant. This is 
in contrast to the collinear antiferromagnets and above considered case J2 > 
where the RG flow can be continued to arbitrary large A. 

To generalize the presented analysis to finite A^ we consider the perturbative 
2+e RG approach (see Appendix^]). For simplicity, we restrict ourselves to the 
RG equations with the space-time isotropy, i.e. J.f = an d use the dimension 
regularization scheme which conserves gauge symmetry. This regularization 
neglects the ultraviolet cutoff originating from the lattice. The structure of the 
flow for N = 3 is shown in the coordinates a = p® n / p° ut — 1, f3 = A/(2n 2 p^ ut ) 
in Fig. 2. The RG trajectories are characterized by the ratio M^ G /p out where 



M RG = — ^ (37) 

M * ~ N\p m /p out -l-a c \("-mN-D V) 

is the generalization of the absolute value of the mass M x of the field in- 
troduced above to finite N. The flow along the trajectories is parametrized by 
A/ p ou t- The RG flow contains a line of trivial fixed points (3 = and two non- 
trivial fixed points. The non-trivial gauge-symmetric fixed point (0(2) GS) 
has coordinates a = —1, f3 = 1/(N — 2) and it is infrared-unstable in both 
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directions. This point belongs to the line a = — 1 where the mass M x van- 
ishes and the theory is gauge invariant. The second point C + has coordinates 
a = a c = (N - 3)/(JV -1),/3=(N- 1)/(N - 2) 2 and it is infrared-unstable 
in one direction (a = a c ) and stable along the other direction (critical line) 
for arbitrary N. Hence, it determines an infrared behaviour near the phase 
transition, corresponding to the bare parameters on the critical line. The cor- 
responding critical exponents cannot be determined reliably by the 2 + e RG 
approach [ID] . Moreover recent investigations have shown the absence of stable 
fixed points in three dimensions at small enough N [HUH]- Nevertheless, not 
too close to the quantum critical point the pseudo-universal behaviour is ex- 
pected, which is correctly described by the 1/iV-expansion. The corresponding 
exponents in the first order in 1/N, (3 = 1/2— 6/(tt 2 N) and v= 1 — 16/(7r 2 iV) 
(see Appendix [D]) coincide with those for the three-dimensional classical phase 
transition determined in the framework of Landau- Ginzburg phenomenology 
[I]. For N = 3 the values of the exponents are in reasonable agreement with 
the pseudo-critical exponents obtained within the non-perturbative RG [6]. 

Two of four regions of the flow diagram of Fig. 2 correspond to the symmetric, 
and two to symmetry broken state. In the regions with a < a c the mass 
M RG determines a crossover from the "massless" (the part of the region near 
0(2) GS fixed point) to "massive" behaviour of the propagators. On the other 
hand, at a > a c where M RG cannot be arbitrary small, it determines a change 
in the RG flow from the universal massive regime at A < 16M RG to the non- 
universal behaviour at A > 16Af RG , the latter essentially depends on a cutoff. 

Similar to the N = oo results, the RG flow at a < a c can be continued to 
arbitrary large A, while, in the right-hand side of the phase diagram (a > a c ) 
a critical value of the scale 



(A*) RG = 2vr [^(+oo) - HPiJPout ~ 1)] M X RG , (38) 



with the function J- {a) defined in Appendix O exists. (A^) RG is proportional 
to the mass M RG defined in Eq. fl37|) . the proportionality factor being non- 
singular. In this case the RG flow cannot be continued to the region A > 
(AJ) RG . As discussed above, for A > (A^) there are no bare parameters 
corresponding to the renormalized Pin iOU t- The value of A is, however, bounded 
by A < A phys where A phys is determined by physical conditions of applicability 
of the continuum model (for frustrated spin systems A p h ys a -1 ). For A p h ys < 
(A^) RG the action QEAF, Eq. ([7]) is applicable in the whole admissible range 
A < Aphys- In this range the QEAF is renormalizable in the limit N = oo, but 
already in the first order in 1/N this property may be lost, cf. next Section. 
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Fig. 2. The renormalization group flow for N = 3 and T = obtained in the per- 
turbative 2 + e RG approach in the coordinates a = pf n / Pout ~~ 1> P = A/(2"7r 2 /9o Ut ). 
The critical (bold) line separates the symmetric (upper part) and symmetry broken 
(lower part) phases; the fixed point C + controls the (pseudo-) critical behaviour 
near the phase transition. The mass scale M^ G vanishes at the long-dashed line; 
0(2) GS is the gauge symmetric fixed point. The bold line a = separates the re- 
gions with trivial (a < 0) and nontrivial (a > 0) ultraviolet behaviour. "Massless" 
and "massive" denote universal regimes with Mf G < A/16 and Mf G > A/16, 
"non-universal" stands for the non-universal regime with M^ G < A/16. Short- 
dashed lines show the crossovers between these regimes. Symbols correspond to the 
bare parameters of TLAF with different spin values (see Sect. 8) 

7 Correlation length at finite temperatures 



In this Section we consider the application of the developed 1/A^-expansion to 
the calculation of correlation length at finite temperatures. 



7. 1 Correlation length in the renormalized classical regime in two dimensions 



In the renormalized classical regime (on the ground-state-ordered side) we 
expect the exponentially large correlation length for T <ti p ou t, similar to the 
collinear case [28] • At N = oo we obtain 



t _i c out 27Tp ou t /on\ 

6v=oo = r = — exp nt . (39) 
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Due to the exponential dependence, leading corrections to the correlation 
length in the first order in 1/N come from momenta and frequencies which 
are much smaller then the scale determined by the temperature. This enables 
one to take into account only terms with zero Matsubara frequency and use 
T/c out as an ultraviolet cutoff for the integrals over momenta. The calculation 
of the logarithmic corrections follows the same lines as in Chubukov et al. |28j . 
The correlation length in the first order in 1/N is given by 

5r 2 1 

r = r + — + Tr ^S(k = iro ,u n = 0), (40) 



where 



1 R A - R x , . f (f'i> ,- <2i 



K -lw? GM - (41) 



The quantities R\ t A are defined by 

Rx = J (03 G o(p)Si(p); Ra = J ^^(p)E 2 (p), (42) 

the self-energies E 12 being defined by Eqs. (126|) and (1271) . In particular, E x 
contains the following logarithmic term 

- 3r 2 T 

Si(k = tr , u n = 0)^ -J 17 -£ lnln . (43) 

2N r c out 

In contrast, S 2 does not contain singularities. Nevertheless both and A 
fields contribute logarithmic terms to the mass correction 5r 2 . Summing up 
all the contributions we obtain 



: J . 
2N 



-RC C °ut / NT 



1 ^_ 1 ln(2p out /p il 



AT iV 2p out /p in - 1 



2vrp 



out 



(44) 



where the last term in the square brackets is due to the field fluctuations. 

In fact, when we turn to smaller values of N, both the exponent term and its 
prefactor in Eq. f)44p are expected to be corrected by higher orders in 1/N. 
Instead of direct calculation of these corrections, it is possible to obtain their 
exact values with the help of RG. For the correlation length we find 

^ « P( Pi J Pout ) C -f exp {2vr [^(p in /p out - 1) - f(a c )} Mf /T} , 

~" I>: ' U,J (45) 
where k = (6N 3 - 27N 2 + 32N - 12)/[2(N - 2) 3 (2N - 3)], F{a) was defined 
in Appendix [Cj P(x) is some function (we do not present here its explicit 
form). We see that the A^-dependence of the correlation length fj4"5]) is quite 
complicated, in comparison to the square lattice case 
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7.2 Correlation length at finite temperatures above the quantum critical point 



From Eq. (ITT)) we obtain that the inverse correlation length above the quantum 
critical point Ji = (Ji) c in the zeroth order in 1/N is linear in temperature, 

r 1 = n, = er, (46) 

where 

6 = 2 In — 1 = 0.962. (47) 

This behaviour of correlation length is modified by the 1/iV-corrections, as 
given by Eq. (1401) where 5r 2 now has an additional term coming from the 
renormalization of the critical value {J\) c 

5 ~2 = Ra-Ra(T = 0)-Rx + R x (T = 0) 

K = a/5/ (87rTG J± ). The general result for the correlation length has the form 

7 (T/A phys ) + W (16M M /A phys ) + U(T/M^ 16M M /A phys 



N 



r^er/Ct i + 

iV 

(49) 

where W and U are some functions, A phys is the upper cutoff scale introduced 
at the end of the Sect. 6. The 7— contribution comes from the fluctuations of 
the A fields (Si and R\ contributions) and coincides up to the factor 3/2 with 
the similar correction to the correlation length of the collinear antiferromagnet 
[28], 7(0) = 0.3560 to first order in 1/N. The function W corresponds to the 
renormalization of the out-of-plane spin- wave velocity by the zero temperature 
fluctuations, c r out = 1 + A^" 1 ^ in the first order in 1/N. This renormalization 
arises because of the absence of the Lorentz invariance leading to the difference 
of the renormalized and bare out-of-plane spin- wave velocities (the latter is put 
to unity). The remaining term U comes from the A field fluctuations (S 2 and 
Ra contributions) and leads to the non-trivial temperature dependence of 
One can expect the universal behaviour of the £ in the regimes M M <C A p h ys /16 
or |M M | ^> A phys /16, where the sum W + U depends on T/M^ only. Due to 
the finite mass of the A field this may lead to a crossover in the temperature 
dependence of £. The latter is in contrast to the LGW-model approach where 
no current terms leading to additional massive field in the quantum critical 
region are considered. 

Consider first the "massless" regime 

M T>X < A phys /16. (50) 

For frustrated lattice systems this implies 16M^ «C a -1 and can be realized 
only in the left part of the phase diagram (i.e. for p in < p out , Xm < Xout)- 
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Condition (!50l) also corresponds to an abstract field theory where the cutoff 
parameter should be chosen larger than any other scale. The asymptotics of the 
integrands allow to conclude that U has a finite limit for A p h ys — > oo. The same 
assertion holds for W. This is in accordance with the renormalizability of the 
theory at finite temperatures. The temperature dependence of the correlation 
length for infinite A p h ys was computed numerically and shown in Fig. 3 for 
M T /M X = 1.095. In a broad temperature region T ^> M T ^ X we obtain quadratic 




0,0-1 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 — . 1 1 

2 4 6 8 10 12 

T/M 

T 

Fig. 3. The temperature dependence of the coefficient C = 0.3560 + W + U 
which determines the correlation length in the quantum critical regime, 
f- 1 = 6T/(£ ut [1 + C/N] for Mp « Aphys/16 (massless regime) and 
M T /M X = 1.095 

behaviour of the inverse correlation length. For T m M TjX the crossover to the 
low-temperature regime 

C\T - 0) = (l + 2^60) (51) 

occurs. Note that in the limit T — >• the contributions W and U cancel each 
other, and the field does not contribute to correlation length. 

In the case of large enough \M T>X \ 3> A phys /16 (massive regime) the theory is 
universal at small enough temperatures T <C A p h ys /16 \M TjX \ in both the 
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right and left part of the phase diagram. Expanding the function U in Eq. (149p 
to first order in l/M TjX we find that it contains a non-universal contribution 
proportional to A phys , 

^ = 7T~2 (t7 n~) Vy- ( 52 ) 



157T 2 \M X M T 

This contribution is, however, cancelled by the similar term in the function W 
which defines the out-of-plane spin-wave velocity renormalization. The result- 
ing universal expression for the correlation length in the first order in l/M XjT 
is 

i / 1 „ ,r , a „„, mnlT / 0.3560 0.0598 T 0.0473 T 

(53) 

In Eq. (1531) we also find both linear and quadratic terms in T. Contrary to the 
case |Mj. )T | <C A p h ys /16, the T 2 -term is small in comparison with T-term. 

For an intermediate region |Ma. jT | ~ A p h ys /16, which corresponds to the crossover 
region from small to large M M (in the left part of the phase diagram) 
or to the non- universal regime (in the right part), we obtain the correlation 
length dependent on A phys . At low temperatures T <C A phys /16 we have £ _1 = 
QT/c r out [1 + (0.3560 + xT)/N] with x being some non-universal constant de- 
pendent on M^/Aphys- At 161-^1 ^> A p h ys the value of x is small and universal, 
see Eq. ( 1531) . In the right side of the phase diagram for 16|M M | ~ A c > A p h ys 
absolute value of x is slightly larger then in the left part and main contri- 
bution to t^~ l (T) is given by the field fluctuations with momentum and 
frequencies of order A phys . 



8 Application to TLAF 



Now we consider application of the obtained results to TLAF. The bare pa- 
rameters corresponding to TLAF can be obtained using asymptotics j\ (A — > 
0) — > Pout/A and a (A — > 0) — > pin/pout — 1 of Eq. (1C.2I) and extending the 
resulting RG trajectories to larger A < A pnys - The fully renormalized spin 
stiffnesses and susceptibilities can be obtained from the spin-wave results [31] 

2 / 0.291\ V3JS 2 ( 0.125 
Xout = - f- T 1 > Pout = 1 



9y/3J V 2S y Fom _ 4 V 25 , rj 

2 / 0.448 \ V3JS 2 ( 0.678^ ° 

Xm = 1 777T ), Pin = — = — 1 



9V3J V 2S J ' rm 2 V 2S 

rearranged according to the scaling laws p in = Am 1 / 2 , p out = Am 1 / 2 (I — 
Bm) (and similar for spin susceptibilities), where A and B are spin-dependent 



19 



constants, and m = 1 — 0.261/S 1 is the sublattice magnetization. The bare 
parameters at A = A phys corresponding to various spin values are marked by 
different symbols in Fig. 2 (the parameter A phys is put to 0.3a -1 ). The open 
triangular, square and circle denote values of spin stiffnesses (right-hand side 
of the phase diagram) and susceptibilities (left-hand side) for S = 1/2, 1, 3/2, 
respectively. Fig. 2 also shows the parameters corresponding to spin values 
S < 1/2 which approach QPT at S = 0.261: the filled triangular, square and 
hexagon correspond to spins S = 0.267, 0.29, and 0.36. The dependence of 
both (A^) RG and M RG on spin is almost linear, (A*) RG ~ 6.365 - 0.34 and 
16M RG ~ 5.075— 0.67; these quantities decrease when approaching QPT. Near 
QPT we obtain (A^) RG ~ 16M RG ~ a~ l . To apply the continuum description 
to this case the inequality A p h ys C a -1 ~ 16M RG must be satisfied, so that 
the result (|53|) for the correlation length is applicable. 



9 Conclusion 

To describe frustrated spin systems we have considered a generalization of 
the previously derived quasiclassical 0(3) x 0(2)/0(2) model to the quantum 
case. The corresponding model includes all the relevant terms, some of which 
vanish in the classical limit. Additional contributions contain both space and 
time components of the 0(2) current interaction term. In the case of the TLAF 
they are proportional to the square of the projection of plaquette spin current 
onto the axis perpendicular to the plaquette. 

To investigate the universality of the finite-temperature properties we have 
constructed the 1/A^-expansion by introducing a generalized N > 2 spin- 
component O(N) x 0(2)/0(N — 2) x 0(2) model. Performing the Hubbard- 
Stratonovich decoupling and introducing Lagrange multipliers for the con- 
straints have led us to the action which contains fluctuating scalar and massive 
vector fields. The former field is responsible for kinematic interaction which is 
below its upper critical dimension. The latter field describes dynamic interac- 
tion which is above its upper critical dimension. 

We have obtained the flow diagram of Fig. 2 and temperature dependence of 
the correlation length of the frustrated spin systems within the 1/A^-expansion 
and 2+e renormalization group approaches, Eqs. (j4"5l) . fl5T|) . and (1531 . The flow 
contains massless, massive and non-universal regimes. In the first two regimes 
physical properties at finite temperature are universal functions of ground- 
state spin stiffnesses, susceptibilities and order parameter. On the other hand, 
in the non-universal regime, the dependence on the regularization procedure 
(e.g. cutoff parameter) cannot be removed by renormalization. When the sys- 
tem is further tuned across the non-universal regime, the continuum model 
becomes inapplicable at scales A > A c , where A c is given by Eq. (I3"S1) . 
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The calculated temperature dependence of the correlation length £ in the 
renormalized classical regime is given by Eqs. (jUJ) and (I4"5l) . In the quan- 
tum critical regime and small enough (massless regime) the correlation 
length £ at very low temperatures shows the dependence £ oc T" 1 , similar 
to the collinear case [28]. At larger temperatures this dependence changes to 
£ oc T~ 2 . Contrary to the concept of the absence of any scales (except temper- 
ature) at the quantum critical point, M M leads to a crossover in temperature 
dependences of physically observable quantities. For large mass (mas- 
sive regime) the temperature dependence of is almost linear with small 
quadratic corrections. These corrections increase and become dependent on the 
cutoff procedure approaching the non-universal regime of the phase diagram. 
As it becomes evident in 1/A-expansion, the (non-) universality is determined 
by the fluctuations of vector field. 

The effect of the vector field is to some extent similar to that of vortices existing 
in fact only for N = 3 [32]. At low temperatures the vortices are bounded and 
therefore not expected to contribute substantially to thermodynamic proper- 
ties. However, at higher temperatures their effect should be considered. The 
detailed investigation of the contribution of vortex configurations and their 
relation to the results obtained in the present paper is the subject of future 
research. Besides that, real physical systems are not truly two-dimensional, 
and effects of interlayer exchange are also to be studied. In the presence of 
this exchange the long-range order is present at finite T < T^ ee i, the calcula- 
tion of the Neel temperature TN ee i should include both spin-wave and vortex 
contributions. 
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Appendices 

A Reduction of the action to the 0(N)/0(N — 2) manifold 

To integrate out rotations of the vector fields e 3 , . . . , e^ for fixed ei and e2, it 
is convenient to write down the action in the matrix form. Let z-th column of 
the orthogonal matrix R(x) be constructed from the components of the vector 
ej(x), where i = 1, . . . , N. Then Eq. (jBj) for q 2 — 0, p = pi = p 2 takes the form 

S gen = \jdx (-p Sp {(R T d T R) 2 } - qi Sp {P(R T VR) 2 }) , (A.l) 
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where P is the projection operator, Pn = P22 = 1 and Pij = in all the other 
cases. When ei, e 2 are fixed, the rotation of the vectors e 3 , . . . , e^v is described 
by the transformation 

R(x) 1 — ► R(x)V T {x), (A.2) 

where 

,h , 

V{x) = , (A.3) 



V'(x) 



so that V T V = 1. 



We introduce the following parametrization for R. For each matrix R we con- 
sider the matrices Rq and V which satisfy the relation R = RqV t . Thereby 
we can factorize the group of rotations O(N) in the factor set 0(N)/0(N — 2) 
and the group 0(N — 2). Then we substitute this representation of R into 
(jA.ip . Note that the contributions coming from V in the terms with spatial 



derivatives vanish since they are cancelled by P. As a result we obtain 

S gen [R] = S gcn [i? ]+| J dx (- Sp {{V T d T Vf} - 2 Sp {(P'R^d T R P'fV T d T V} 

(A.4) 

where P' = 1 — P is the projector complementary to P. It is important that the 
interaction between the fields Rq and V is equivalent to the quadratic interac- 
tion between their currents Rq8 t Ro and V T d T V. Besides that, the kinematic 
interaction which originates from the integration measure is also possible. 
However, as we will show below, all the contributions from the measure of in- 
tegration are trivial. Now we replace the variables V{r) by V T d T V . Since the 
space of antisymmetric matrices is flat, the integral over V T d T V is gaussian. 
Calculating this integral we arrive at the effective action 

S Q eaf[R ] = S gcn [R ] + ^JdxSp {(P'R^d T RoP') 2 }, (A.5) 

with the order parameter defined on the 0(N)/0(N — 2) manifold. In Eq. 
(IA.5I) we return to the e a representation and obtain the action (J7J) of the main 
text with J{ = 2p, Jf = Jf = q x , J T 2 = p, Jf = Jf = 0. 

Now we show that the measure of integration gives a constant contribution 
to the action (J7|). The transformation of the invariant measure dfi(R) of the 
rotation group has the form 

dfi(R) = dfi(RoV T ) = f(Ro,V)d/2(V)du(R ), 

where dv is the left-invariant measure on the 0(N)/0(N — 2) manifold and 
f(Ro, V) is the corresponding Jacobian. From the invariance of the left side of 
this identity with respect to the transformation R \ — > gRh we have f(Ro, V) = 
f(gRo, h T V). This means that the Jacobian / is constant and the kinematical 
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interaction between Rq and V is absent: 

dfi(R) = Const dn(V)dv(Ro) 

It remains to show that the Jacobian of the transformation from V(r) to 
V T d T V is constant. To this end we write 

Y[dfi{V{r)) = l[H{V T d T V)d{V T d T V). 

T T 

The invariance of the left side with respect to the gauge transformation V i — > 
h(r)V leads to H(V T d T V) =const. 



B Global stability of the saddle point 

In this Appendix we show that the saddle point defined by Eq. (ITT)) gives 
main contribution to the partition function (II 2p for N — > oo. Actually we 
want to use the following mathematical theorem for the saddle point method 
[l3|H]: if the minimum of the function Re {S^ff} along the integration contour 
is achieved only at the saddle point P, then in the limit N —>■ oo the integral 
( fT2l) equals to the contribution from this point. We suppose that the system 
is placed into a finite volume with periodic boundary conditions. In this case 
we can safely integrate out the gaussian fluctuations of a a in ( fl2l) . The result 
of integration equals Eq. (fT2|) apart from a a = 0. We assume also that the 
quantity — > is nonnegative since it equals to half of the bare in-plane 
spin stiffness (for fi = x) or susceptibility (for fi = r) which are supposed to 
be positive. 

Consider first the case J% > where the field is real. We can rewrite the 
real part of the effective action (lT3"j) as follows: 

2Re^ eff [A,A] =Re^\n(l + iK[A]- 1 L) + splnK[A], (B.l) 

where K[A] = (Jf - J%){id^) 2 + J 2 M (i9 M - o~ v A^) 2 + v 2 is a positively defined 
operator, A^ and L a/3 (x) = — i(\ a p(x) — v 2 5 a p) are arbitrary real functions 
of space and time variables. The operators i^[v4] _1 and L do not commute 
in general case. However, since all the eigenvalues E n of the product of the 
positive and hermitian operators i^[A] _1 L are real, we obtain for the first term 
in Eq. flBTTT) 

Re sp In (l + iKIA^L) = - £ Re In (1 + E 2 n ) ^ 0. (B.2) 

It can be shown that the sign of equality in this condition stands only for 
L = 0. 
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The second term in ( IB.lj) is expressed as a functional integral over complex 
scalar field f(x) 



(B.3) 



spin K[A] = -2 In (J D[<p* , <p] exp |- J dx(p*(x) [(Jf - J%)(id t 
+J^ l d^-A,f + v 2 ] <p(x)}). 

Using transformation of variables <p(x) = (p r exp(ix) we obtain 



splnf^fA] = —2 In (J (p r D[(p r ]G[A, cp r ] exp |— J dx(p r (x) (—Jid 2 + v 2 ^) (p r {x) 

(B.4) 

where 



G[A, cp r ] = J D[ X ] exp {- J dx<p r (x) 2 [(Jf - 4)(d, X f + 4(A, - d.xf] } 

(B.5) 

For a fixed <p r the maximum of the functional G[A, (p r ] is reached at A^ = 0. 
Therefore, it follows from Eq. flB.4j) . that the minimum of sp In K [A] is achieved 
only at A^ = 0. Consequently, using Eqs. (|B.2[) . (IB. 1|) and aforementioned 
theorem we conclude the global stability of the saddle point (A, A) = (v 2 , 0) 
defined in (TT7I) for > 0. Even if there are some additional saddle points, 
perhaps with coordinate dependent A^ix) and A(x), their contributions to 
(Tl2l) are exponentially small at N — > oo. 



For J% < 0, when the field A^ = iB^ is imaginary, we have 



2Re5' efr [A,A] = Respln(l + iEr 1 y) +splnC7, (B.6) 



where U = Jf (z^) 2 + 1 J^B^ + v 2 , V = L+| J|V(^ S m + B^). The above 
considerations lead us to the stability of the saddle point (A, A) = {v 2 , 0) in 
this case too. 



Note that in the same manner one can prove the global stability of saddle 
points which arise in the 1/iV-expansions of some nonlinear models, e.g., vector 
O(N) NLcr model and CP^-model. Main requirement for the applicability of 
the method is the stability of the action after removing constraints. 
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C Perturbative 2 + e RG approach 



The one loop 2 + e RG equations for QEAF in the case J* = J% (a = 1, 2) 
and T = where obtained in Ref. [5] 

dx a + 5 
-ex + N — , 



din A z /q j\ 

di^A = — {N - 1){a - ac) > 

where x = p° out /(K 2+e A e ) >0,a = > -1, A D = [2 D_1 7r D / 2 r(D/2)] _1 , 

a c = (N — 3)/(N — 1), bare parameters p^, p°ut are related to Jf 2 through 
Eq. (0). The solution to Eqs. (jCl) is 

x = -i ^- jfq, 

x 1 + a n (c.2) 

A e /C = / + — 
n 

where X is the RG invariant (X > for the symmetric phase and X < for 
the symmetry broken phase), /C > is an integration constant, 

v /2N-l\ /N-2 

^( Q) = W i + Qc )]-^r(— )r(— 



N 1 2A-1 a, 

[ V (l + a)} *-i 2 Fi 



A- 1' A- 1' A- 1 ' 



(C.3) 



where 2-P1 is the hypergeometric function, 77 = 1 for a > a c and —1 for a < a c . 
A dimensionfull constant K, playing the role of the crossover scale of the RG 
flow can be expressed through a renormalized parameters of the ordered state 
Pin,out for X < using the relations x = p out / (K 2+e A £ ), a = p h J p out - 1 in the 
limit A -> 0. For e = 1 we obtain K, = l/(7r 2 AM RG ), M rg being defined in 
Eq. d37D. 



An important property of T{a) is that it has a finite limit for a —>■ 00. Hence, 
the RG flow cannot be extended to the region A > (A^) RG , 

(A * )RG .p±i^M) 1/e . (a4) 

Note that this restriction of the RG flow should not be confused with pecu- 
liarities of the RG flow originating from the inapplicability of perturbation 
theory in the symmetric phases. Indeed, the scale (A*) exists even in sym- 
metry broken phase. For X < the scale (A^) RG is expressed through the 
renormalized parameters of the ordered state, as in Eq. 
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D Critical exponents in dimension 2 < D < 4 near QPT 



In this Appendix we consider the calculation of critical exponents of QPT in 
space-time dimensions 2 < D < 4 under the assumption that the transition is 
of second order. 

First we determine the critical exponent for the order parameter a = |mi| = 
|m.2 1 in the symmetry broken phase. The expression for a can be obtained 
from the identity 

N r d 3 v 

Mx)e p (x)) = £ J j^GUp) = <W (D.1) 
which follows from the constraint e^^ = 8 a p. Using Eq. ([22]) we obtain 



a 2 



1- / -^G (p)+Rx-Ra 



(27T) 



;i + f a -Fa), (d.2) 



where 

3 f d D q h(q) 1 f d D q 



1 /■ ^ 3G S (p). „ _ 1 M V Gj^Jf^fi-i^, (D .4, 



iVJ (2tt) 3 n(p) A^J (2tt) 3 
and 

K D S D (3 - D) D _ 6 K D S D 2 D _ 6 

J H<?) = 73 9 ' J 2 (^) = — 73— (l V - > ( D - 5 ) 

T 2 (D/2) 

bD ~ r(D - 1) sin [7r(D - 2)/2] ' ( j 

In the critical region Eq. ( 1D.2I) takes the form 

2 ^ A Jf \ 13(5 -2D), / K D S D A D - 2 \\ 

where Jf is the critical value of the coupling constant J% which depends on 
J^. It follows from flD.71) that the critical exponent /3 for the magnetization in 
the first order in 1/N is given by 

pi»( D ) = I + I JL^D (D.8) 
2 N 2Sd D — 2 K 1 

The critical exponent v can be determined from the scaling behaviour of the 
spin stiffness 

p out oc(l- JrM) (D - 2)i/ , (D.9) 
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which is determined by Eq. ( !30|) with o"o being replaced by renormalized order 
parameter a. For J x - J{ < |M M | we obtain from Eqs. (1221 . (I24"|) . and ([25]) 

(D.10) 

As a result we obtain 
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